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ABSTRACT 

Galaxies become red and dead when the central supermassive black hole (SMBH) 
becomes massive enough to drive an outflow beyond the virial radius of the halo. We 
show that this final SMBH mass is larger than the final SMBH mass in the bulge 
of a spiral galaxy by up to an order of magnitude. The M — a relations in the two 
galaxy types are almost parallel (M oc er 4+ ^, with (3 < 1) but offset in normalization, 
with the extra SMBH mass supplied by the major merger transforming the galaxy 
into an elliptical, or by mass gain in a galaxy cluster. This agrees with recent findings 
that SMBH in two Brightest Cluster Galaxies are ~ 10 x the expected M — a mass. 
We show that these results do not strongly depend on the assumed profile of the dark 
matter halo, so analytic estimates found for an isothermal potential are approximately 
valid in all realistic cases. 

Our results imply that there are in practice actually three M — a relations, cor- 
responding to spiral galaxies with evolved bulges, field elliptical galaxies and cluster 
centre elliptical galaxies. A fourth relation, corresponding to cluster spiral galaxies, 
is also possible, but such galaxies are expected to be rare. All these relations have 
the form ATbh = C n a^ , with only slight difference in slope between field and cluster 
galaxies, but with slightly different coefficients C n . Conflating data from galaxies of 
different types and fitting a single relation to them tends to produce a higher power 
of a. 
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1 INTRODUCTION 

In the past decade, observations of large numbers of galax- 
ies have revealed important correlations between the prop- 
erties of galaxy spheroids and the supermassive black holes 
(SMBHs) they host. Two of these relations are particularly 
interesting: the black hole mass - host spheroid velocity 
dispersion (M — a) relation (|Ferrarese et al.ll2006l) and the 
black hole - bulge mass (M — Mb) relation ( Haring fe Rbd 
120041 '). Together, they strongly imply that SMBHs coevolve 
with their hosts, and probably affect each other's evolution 
through some form of feedback. 

There have been numer ous attempts to explain these 
relations, both analytically (fsTlk fc Reesl [l99g| ; ikind 12001 
120051 : iMurrav et al.l 120051: [Power et all 1201 lT) and numer - 



Tombesi et al"1l2010ll?h a nd outflows^ (e.g. lRupke fe Veilleuxl 
2011 : ISturm et alT 201 lh . In this picture, the wind and out- 



ically (e.g. |Pi Matteo et all 120051 : iBooth fc Schavd I20u9h 



The AGN wind feedb ack model (|Kind 12001 120051 . l201ol : 



I Zubovas fc Kind" |2012| . and references therein) is particu- 
larly promising, as it explains not only the observed cor- 
relations, b ut also other obser v able properties of galactic 
winds (e.g. iPounds et al.ll2003l l?l: IPounds fc Vaughanll201ll : 
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flow shocks initially occur fairly close to the SMBH, and 
are efficiently Compton-cooled. The outflow retains only the 
ram pressure o± L^dd/c of the black hole wind and has ki- 
netic energy lower by a factor a/c ~ 10 -3 . Only when the 
SMBH mass approaches the critical mass M CT ~ 3.7x10 s M 
(see equation below) do the shocks move further away. The 
resulting geometrical dilution of the quasar radiation field 
now makes Compton cooling inefficient. The shocks are no 
longer cooled, and the ou tflows become energy-driven. 

In two recent papers (|King et al.ll201ll : IZubovas fc Kind 

I2012T ) we investigated the properties of large (kiloparsec)- 
scale outflows in galaxies. These must be energy-driven, 
i.e. the outflow kinetic energy rate E ut = M out «out/2 — 



1 For clarity, we distinguish between the black hole wind, the gas 
coming directly from the accretion disc around the SMBH, and 
the outflow, the outward movement of swept-up interstellar gas 
from the host galaxy. These two components undergo respectively 
reverse and forward shocks on each side of the contact disconti - 
nuity separating them. See the Figure in lZubovas fc Kind i2012T ) 
for more details 
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0.051/Edd, where I/Edd is the Eddington luminosity of the 
driving quasar. We showed that energy-driven outflows are 
able to clear galaxies of gas, turning them into red-and-dead 
spheroids. 

However, in these papers, we only briefly touched upon 
the question of how much the SMBH mass grows as it expels 
the outflow from the galaxy, an d this remains unce rtain. 
There are two discrepant claims: I Power et al.l l|201lf ) used 
an energy argument to conclude that the SMBH only needs 
to grow its mass by ~ 40% t o expel the surrou nding gas 
past the virial radius, while in iKing et al.l (|201ll ) we found 
that an outburst duration of ~ 10 s yr is required, which 
would lead to the SMBH mass growing by about an order 
of magnitude. 

We return to this problem in this paper. We also con- 
sider the difference between evolved spiral and elliptical 
galaxies in more general terms, as well as the dependence of 
the outflow properties on the galaxy environment. We find 
only a small uncertainty in our results stemming from the 
assumption of an idealized background potential (Section[2]). 
By considering the distribution of energy in various compo- 
nents of the outflow we resolve the discrepancy between the 
two claims (Section |3J. We conclude that SMBHs in ellipti- 
cal galaxies should grow to a mass a few times higher than 
in spirals for a given value of a. We consider the effects of 
gas depletion and replenishment in a galaxy and find that 
galaxies in cluster environments should have slightly more 
massive SMBHs than their counterparts in the field (Section 
[4|. We summarize and discuss our results in Section [5] 



2 OUTFLOWS IN A GENERAL POTENTIAL 

In previous work on the wind feedback model, we adopted 
an isothermal background potential, which has the use- 
ful property that the weight of the gas is independent of 
radius. To check the validity of this assumption we now 
consider what happens to an energy-driven (large-scale) 
flow in a more reali stic halo. As an e xample we con- 
sider an NFW halo jNayarro et all I1997T ). We note that 
iMcQuillin fc McLaughlin! (|2012T ) have considered this prob- 
lem assuming momentum-driven flows for all radii (also for 
Hernquist and Dehnen-McLaughlin potentials) and shown 
that the differences from the simple isothermal treatment 
are fairly small. 

One particular difference between isothermal and NFW 
halos that is important for us is the varying velocity disper- 
sion. In general, all non-isothermal halos have varying a, 
usually with a increasing from the centre out to some ra- 
dius r pea k and decreasing beyond it. As the weight of the gas 
swept up by the outflow depends on the value of a through 
the equation governing its density, it is no longer indepen- 
dent of radius, as is the case for an isothermal distribution. 
Specifically, the weight of the gas increases as the outflow 
moves from low radii toward r poa k and decreases afterwards. 
This means that an SMBH is able to inflate a rather large 
(but still with R < i? pca k) outflow bubble even though its 
mass may still be too small to completely drive the gas out of 
a galaxy. On the other hand, once the outflow passes i? pca k, 
it becomes easier to expel the gas. 

In the present calculations, we consider an NFW halo 
with a concentration parameter |Navarro et alj 119971 1 c = 



10, scale radius i? s = 25 kpc (so that Ry = 250 kpc) and 
virial mass Mv = 3.3 x 1O 12 M . This gives <7 poa k - 200 km 
s _1 at 7? pcak ~ 2.167? s ~ 54 kpc. We take the initial SMBH 
mass to be M = 3.68 x 1O 8 M , eq ual tq the c ritical SMBH 
mass found in previous papers (e.g. iKinakoiCf ). We compare 
this NFW halo with an isothermal halo with a = 200 km/s. 
We do not truncate either halo beyond Rv, but allow them 
to continue to infinity. 

The equation of motion for an en ergy d riven -shell in 
an isothermal potential was derived by iKind (|2005l ) and its 
solutions investigated in IKing et al l (|201ll ). Here we derive 
a generic equation of motion, valid for any spherically sym- 
metric mass profile, provided that the gas density traces the 
background potential density, i.e. M g (R) = f g M (R) for any 
R, where / g is called the gas frac tion. We start with equa- 
tions (15) and (16) of lKind (|2005l ): 
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Using eq. ^ to eliminate P from eq. ((2| and expanding 
the terms with derivatives leads to a generalized equation of 
motion (we have dropped the (R) identifier from each mass 
for compactness): 
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Here, M = dM/dR x R and M = dM/dR x R + R? 



d 2 M/dR 2 . This equation reduces to eq. (18) of lKind (120051 ) 
in the case of an isothermal potential. 

For an NFW potential the equation cannot be solved 
analytically, so we integrate it numerically. This also allows 
us to follow the growth of the central black hole mass (which 
slightly affects I/Edd on the left hand side of @). Figure 
[T] shows the propagation of four outflows, with f s = 3 x 
10 -3 , 0.01, 0.03 and 0.16 (lines from top to bottom in the 
top panel, respectively), in an isothermal potential (left) and 
an NFW potential described above (right). In both cases 
the SMBH is radiating and growing at the Eddington rate 
for t — 90 Myr = 2tsai • Note that our results are valid 
only where the outflow velocity is i> out < u w — 0.1c. This 
condition is certainly satisfied when the gas density is high. 

We see immediately that in an NFW halo, the outflows 
reach larger radii than in the isothermal case. This is ex- 
pected, as the outflows are pushed past the radius of high- 
est a and so are able to accelerate to higher velocities before 
the central AGN switches off. The outflow velocity is almost 
constant during driving in the isothermal case, growing only 
in response to the growth of the SMBH (v tx M^Q) and de- 
cays approximately as v oc i -1 / 2 afterwards. In the NFW 
case, the velocity is not constant, but the minimum velocity 
of the densest outflow is reached at R ~ -R pea k and is similar 
to w e , as might be expected. 

From the results of the numerical integration it is clear 
that the isothermal potential probably represents a worst 
case scenario for outflow escape. In other words, the required 
AGN activity time and the final SMBH mass obtained using 
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Figure 1. Propagation of energy-driven outflows in isothermal (left) and NFW (right) halos. The four outflows have / g = 3 X 
10 -3 , 0.01, 0.03 and 0.16 for the black solid, blue dashed, green double-dashed a nd red triple-dash ed curves, respectively. The SMBH 
is active for 90 Myr (dashed vertical line in top and middle panels). As shown in lKing et al.l d201lfl . outflow properties do not strongly 
depend on initial conditions (radius and velocity). Top panels: radius of the contact discontinuity as function of time. Middle panels: 
outflow velocity as function of time. Bottom panels: outflow velocity as function of radius. 



the isothermal assumption are upper limits. This is reason- 
able, because the weight of the gas, which is oc a 4 , is con- 
stant for the isothermal halo, whereas it varies and is always 
smaller for the NFW halo, which has a < a pca k = o"i so thcrm. 

The differences between the isothermal and NFW mod- 
els are rather small for high gas density. For lower gas den- 
sities (/ g <C 0.16), the outflow in an NFW potential reaches 
much higher velocities and escapes past Rv much faster. 



However, since the higher densi ty case is re levant for estab- 
lishing the final SMBH mass dKind [20031 ). we expect the 
uncertainty due to the shape of the background potential 
to be modest. This is reassuri ng in that earlier treatment s 
of energy-driven outflows (e.g. lKindl2005l : iKing et alJkoill ) 
remain valid. However the question of the final SMBH mass 
remains. In the rest of this paper, we analyze this question 
in more detail. We use the singular isothermal sphere pro- 
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file of both gas and the background potential. We motivate 
our choice by the fact that analytical calculations are much 
more difficult, and often impossible, for non-isothermal po- 
tentials. In addition, using a different profile does not nec- 
essarily make our results more realistic: the central parts of 
galaxies, which are baryon do minated, are kno wn to have 
profiles similar to isothermal l|Naab et al.l 12007), and even 
their dark matter profiles seem to differ from any known 
analytic profiles. 



velocity; while theoreti cally we expect it to a l ways be similar 
to 0. 1c, observations jTombesi et al.l [2010l F^; IPounds et al.l 
l2003l l?l) suggest that it may vary from ~ 0.03c to ~ 0.2c, 
leading to variations in Rc between ~ 50 pc and ~ 2 kpc 
for a 1O 8 M black hole in a a = 200 km s 1 potential. 

We see that close to the SMBH, outflows are 
momentum-driven. They only become energy-driven once 
the SMBH reaches some critical mas s which allo ws a driven 
outflow formally to reach any radius (King 2005) . This mass 



3 PROPERTIES OF LARGE SCALE 
OUTFLOWS 

The radiation pressure of an active SMBH drives a quasi- 
spher ical wind from the outskirts of its accretion disc |Kind 
|2010| ). which then hits and shocks against the ambient in- 
terstellar medium (ISM). The wind is mildly relativistic 
(» w ~ 77c ~ 0.1c, where r\ is the radiative accretion effi- 
ciency), and so its shock temperature is a few times 10 11 K, 
much greater than the Compton equilibrium temperature 
(Tc ~ 2 x 10 7 K). As we remarked above, the shocked 
wind therefore cools via the inverse Compt on process against 
the p hotons of the AGN radiation field l)Ciotti fc Ostriken 
Il997l ). If the interaction between the wind and th e ISM hap- 
pens within a critical cooling radius (|Kind|20 03[), the wind 
loses most of its original energy and only communicates its 
ram pressure to the ISM, creating a momentum-driven out- 
flow (|Kindl2010h . Outside the cooling radius, on the other 
hand, the wind cannot cool efficiently and most of its energy 
rate _E W ~ r]L^^/2 is communicated to the ISM, giving an 
energy-driven outflow (|Kiriell2005l ; IZubovas fc Kinal2012l ). 



3.1 Momentum-driven flows and M — a relation 

We estimate the cool ing radius Rc by comparing the wind 
cooling timescale (cf. iKin j|2003l ): 
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where M 8 is the SMBH mass Mbh in units of 10 s M , with 
the ISM flow timescale 
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where CT200 is the velocity dispersion in units of 200 km/s, 
«out,m is th e momentum -driven outflow velocity (see equa- 
tion (14) in King 20031 ). and / c = 0.16 is the cosmological 
value of the baryon-to-dark-matter density fraction. Equat- 
ing the two we find 
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Here Vq.i is the wind velocity in units of 0.1c. If we had used 
an expressi on for the energy -driven outflow velocity (e.g. 
eq. 13 from iKing et al.l 1201 ll ) rather than the momentum- 
driven one, the cooling radius would be smaller, so it is the 
momentum-driven velocity that is relevant here. We note 
that the value of Rc depends strongly on the assumed wind 
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This value is very cl ose to the observed M — a relation 
llFerrarese et al.l l2006h . provided that / g ~ / c . Note that 
m iKind (| 2003l). the expression for Al a had an extra factor of 



1/2 in it, because that paper simply assumed an esca pe ve- 
locity ~ a, rather than the full gravitational potential. IKind 
(2005) used a f ull isotherm a l pote ntial, which gives the ex- 
pression 0. In|Kjng_ et al.l (|201ll ). we mistakenly gave the 
incorrect numerical value for M a , although we used the cor- 
rect analytical expression. 



3.2 Clearing the bulge of a spiral galaxy 

Once the black hole mass reaches this value, the outflow can 
propagate to large scales as the SMBH continues to grow. 
The timescale for the outflow to reach the cooling radius is 
~ iflow (eq.EJ), i- e - 2 — 3 Myr for the depleted bulge. Then the 
outflow becomes energy driven and rapidly attains a much 
higher constant velocity 



3/c ' 
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km s 



(8) 



where f' s is the gas fraction in the bulge during the outflow, 
typically of the same order as f c for young galaxies. This 
outflow sweeps the bulge in an additional time 



Rh - Rc 



(i?b,kpc - 0.5) Myr. 



(9) 



For typical bulge sizes of 1 — 2 kpc, this amounts to an- 
other ~ 1 Myr of SMBH activity. Gas expelled from the 
bulge does not necessarily escape the galaxy altogether, but 
it may settle on to t he galaxy disc, which is not affected by 
the outflow direct ly dNavakshin et ail [2012). as part of the 
galactic fountain IjBregmanll 19801 ). 

Overall, this expulsion means that the SMBH must be 
active for ~ 5 Myr after it reaches the M — a mass (JT} , even 
when the outflow propagates in an isothermal potential; in 
an NFW potential, the outflow would propagate faster and 
so the bulge would be cleared sooner. Even accreting at the 
Eddington rate (this may be possible contemporaneously 
with an outflow due to the presence of a massive accre- 
tion disc around the SMBH, which cannot be lifted by its 
wind due to large weight; see Zubovas & Navak shm][2012l ; 
IZubovas et"ai]|201ll ; iNavakshin fc Power! |2010| ). the SMBH 
only grows its mass by ~ 10% during this time. Therefore, 
the value M a (eq. [7} is likely to be the approximate upper 
limit to the mass of supermassive black holes in the bulges 
of spiral galaxies. 
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3.3 Elliptical galaxies and the transition to 
red-and-dead 

Following a major merger, a galaxy changes from spiral 
to elliptical morphology. This means that the size of the 
spheroid grows from a bulge several kpc across, to the 
scale of the whole galaxy, which may have a radius of 
tens of kiloparsecs. Major mergers tend to cause bursts of 
star formation fo llowed by a burst of quasar activity some 
~ 300 Myr later ilSanders et al.lll988l;lGao fc Solomonlll99g|; 
ICanalizo fc Stocktodl200ll ; ICombesll20011 ). This quasar out- 
burst is presumabl y what drives most of t he remaining gas 
out of the galaxy (|Zubovas fc Kinall2012l ). transforming it 
into a red-and-dead elliptical. Noting that the virial radius 
-Rv ^ Rc, we may calculate the duration of quasar activ- 
ity necessary to eventu ally drive the ga s outside the virial 
radius of the galaxy (cf . iKing et alj|2011f ) . Unlike in the sim- 
pler case of a bulge (eq.|§} , we have to consider the expansion 
of the outflow bubble even after the quasar has switched off. 
The result, in that case, is 



7H \v, 



2/3 



yr ~ 2tsai, 
(10) 

where H is the Hubble parameter and tsai — 45 Myr is 
the Salpeter time. We also allow some variation of the gas 
density by including a dependence on / g . If the galaxy is 
almost devoid of gas already, the black hole only needs to 
be active for a very short time to drive an outflow past Rv. 
This is also evident from the top panels of Figure [1] For a 
high gas density f s ~ f c , as may be expected after a major 
gas-rich merger, we see that the final mass of the black hole 
is 



M fln ~ M a x exp(2) ~ 7.5M CT 



(11) 



So the final mass of the SMBH in the centre of a red- 
and-dead elliptical galaxy is expected to be almost an order 
of magnitude greater than that inside a spiral galaxy with 
a red bulge exhibiting low star formation rate and low gas 

density. 

This result disagrees with the estimate in lPower et al.l 

(2011), where we found, using an energy argument, that the 
SMBH only needs to grow for ~ OAtsui in order to impart 
enough energy to the gas to let it escape past the virial 
radius. There are two complications with this calculation 
which give rise to the discrepancy. The first, minor, issue 
is that we used the formal escape velocity from an isother- 
mal potential, w eS c = 2a, to obtain the required energy in- 
put. However, the escape velocity is not well defined for an 
isothermal potential, and the adopted value 2a only moves 
gas to a position with radius a factor 20 larger. As a re- 
sult, gas initially within Ry/20 ~ 20 kpc needs more energy 
input to escape outside the virial radius. Nevertheless, this 
effect by itself only produces a small correction. 

The second issue is more important. As an escape con- 
dition, we compare the virial radius of the galaxy with the 
radius of the contact discontinuity between the wind and 
the outflowing ISM. Even while the quasar is active and 
the outflow proceeds with velocity v e , two-thirds of the in- 
put e nergy is transferred to the ISM outside the disconti- 
nuity (jZubovas fc Kindl2012h . After the quasar switches off, 
this fraction increases above 2/3. By the time the outflow 
stalls, numerical simulations |Zubovas fc Navakshinl l2012l l 



show that the energy content of the shocked wind is only a 
few percent of the total input energy. This happens because 
most of the shocked wind energy is thermal, and the wind 
cools as it expands, transferring energy to the outflow; by the 
time the outflow stalls, the energy per unit mass is roughly 
the same in both the wind and the ISM, so the ISM, which 
contains most of the mass, also has most of the energy. 

All this means that in order to get the inner surface of 
the swept -up host ISM (the contact discontinuity) to expand 
past Ry, the quasar has to be active for at least a few times 
longer, and perhaps an order o f magnitude longer, than es- 
timated by IPower et al.l l|201ll ). This correctio n brings the 
quasa r duration required by energy arguments IPower et al.l 
1)20111 ) into agreement with that required by the bubble dy- 
namics, as here. 

We can now combine the final SMBH mass derived 
above (eq. Il ip an d the bulge stellar mass estimate from 
IPower et all l|201lf ). 



M b < 0.6CM g , vir CT2oo, 



(12) 



where ( ~ 1 is a factor encompassing our uncertainty regard- 
ing the amount of feedback that self-regulated star formation 
produces when compared with the typical gas momentum in 
the bulge, and M g , v i r is the total gas mass inside the virial 
radius of the galaxy. Together with the formulae for mass of 
an isothermal s phere and the virial radius (equations 3 & 18 
in IKing et al.l l|201ll 'l. respectively), they give us the black 
hole - bulge mass relation: 



Mfi n 8 4 

— — ~ 7.5 x 3.7 x 10 a 200 
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0.6C x 2f c a 3 a 20 o 



3 x 10~ 3 C _1 . 

(13) 

This is consisten t with the observed relation (e.g. 
lHaring fc R ix 2004). For spiral galaxies, the black hole mass 
is approximately an order of magnitude lower, but the 
spheroidal component is also smaller by a similar factor, 
so the relation holds roughly in this case as well. 



4 GALAXY ENVIRONMENT AND ITS 
INFLUENCE ON THE M - a RELATION 

So far, in our calculations we have assumed that f g ~ f c 
throughout the evolution of a galaxy. This may be a reason- 
able assumption on scales larger than _Rv, but is unlikely 
to be correct on the smaller scales that most immediately 
affect the properties of the black hole. Galaxy bulges are 
baryon-dominated and so by the time the SMBH begins 
driving large-scale outflows, the fraction of gas remaining 
in the bulge depends mainly on depletion by star formation 
and refilling via mergers or accretion of cold gas from the 
halo. 



4.1 Gas depletion 

When a galaxy forms, its central regions have / g ~ 1, as 
most of the baryons have not yet had time to form stars. 
Subsequently, star formation depletes the gas on a timescale 
tsF- As a rough approximation, we assume that star forma- 
tion persists throughout the lifetime of the galax y, with 2% 
of gas converted into stars on a dynamical time 



y, witn Z7o 
(|Kennicuttl 
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Il998l ). With a scale radius for the bulge of 1 kpc, typical for 
a spiral galaxy, this gives a star formation timescale 

t SF = ^ ~ 2.5 x lO^kpcO-^eo.M yr ~ 5.5t S ai- (14) 

If the SMBH grows from a seed mass of ~ IOMq, it takes 
~ 16£sai to reach 1O 8 M0, by which time the gas fraction is 
depleted to 



' exp 



16t Sa 

tSF 



0.05 ~ 0.35/c 



(15) 



Substituting f' s instead of f B = f c into equation [JJ we find 
that in spiral galaxies we may expect SMBH masses to be a 
factor ~ 2.5 lower than predicted above. 

The estimate above, however, is highly uncertain. At 
least three major effects may drive it to either higher or 
lower estimates for the SMBH mass. The estimate may 
be decreased if the SMBH growth is intermittent, i.e. the 
duty cycle is low, while star formation efficiency is approx- 
imately constant throughout the lifetime of the galaxy. In 
that case, the gas fraction in the bulge may be depleted 
to even lower values than f' g , leading to a smaller critical 
SMB H mass. If, for example , the black hole duty cycle is 
10% (|Schawinski et al J 120071) . then the critical black hole 
mass drops to a very low value of 10 4 Mq. The fact that 
such small black holes are not observed suggests that star 
formation and SMBH growth are temporally correlated in 
galaxies. 

On the other hand, there have been suggestions that 
SMBHs may grow fro m seeds as massive as 10 s M Q (e.g. 
iBegelman et all I2006T ). In this case, the SMBH growth 
timescale is shorter than the 16£sai used above, and so the 
gas is only depleted to f s ~ 0.4, requiring a larger final 
SMBH mass. 

Finally, a third complication may affect the slope of 
the predicted Mbh — o relation. Galaxies with larger a have 
large r effective radii and larger bulges fe.g. lMarconi fc Hunt! 
l200&t ) and so presumably are less efficient in forming stars. 
As a result, the critical value of / g may correlate positively 
with <r, increasing the slope of the mass - velocity dispersion 
relation. This may be important when comparing spiral and 
elliptical galaxies: not only do elliptical galaxies need to grow 
their SMBHs for longer in order to expel the gas past the 
virial radius (see Section |3.3[) . they must do so while the 
galaxy has a higher gas fraction than a spiral galaxy would. 
However, we cannot quantify this effect in more detail with- 
out recourse to a much more detailed model of long-term 
galaxy evolution, which is beyond the scope of this paper. 

A conclusion of this section is that it is likely, but by 
no means certain, that isolated spiral galaxies could have 
central SMBHs that are a factor ~ 2.5 less massive than 
predicted previously. Isolated elliptical galaxies may be sim- 
ilarly offset from the relation given in eq. In addition, 
the slope in the M — a relation for these galaxies may be 
slightly steeper than M oc a 4 due to a correlation between 
the velocity dispersion and size of the galaxy spheroidal com- 
ponents. 
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Figure 2. The four M — cr relations (solid lines) and their com- 
bined effect on observational fits (dashed line) . All solid lines have 
slopes M oc <t 4 and the dashed line has M oc <r 6 . The grey area is 
the ap proximate locus of data points in Fig. 3 of lMcConnell et al.l 

feoiil) . 



inside a cluster, r am pressure and tid al forces strip most of 
its gas away (e.g. [Takeda et aT1ll984M . This process lowers 
the average gas fraction in the galaxy; however, since gas 
is predominantly s tripped from the outskirts of the galaxy 
l|Abadi et alJll999h and star formation is suppressed in lin e 
with gas depletion (|Lewis et al.ll2002j ; iBalogh et al.l 120001 "). 
the limiting central SMBH mass is not significantly affected 
by this process. 

The cluster environment may also provide a replenish- 
ment of gas for galaxies that do not experience significant 
stripping, such as Brightest Cluster Galaxies (BCGs). These 
galaxies sit close to the centre of the cluster potential and 
thus are embedded in a massive hot h alo. Cooling fl ows from 
the halo replenish gas in the galaxy ((Fabian 1994|), sustain- 
ing star formation and feeding the SMBH. If this replenish- 
ment happens on a timescale which is not too different from 
the cluster dynamical timescale td,a = R c \/o c \, then the 
two effects of star formation and cooling flow replenishment 
balance each other at a gas fraction 



<SF 



td,c\ + tSF 



(16) 



Using the same values as above for star formation and R c \ 
1 Mpc and a c \ = 1000 km/s, we find 



/ cq = 0.2 ~ / c 



(17) 



We see that BCGs, and perhaps some other cluster galaxies 
which are moving slowly with respect to the cluster halo 
gas, should have black hole masses governed by the M a of 
equations [TTj and [7] The latter equation may apply to the 
strongly perturbed spiral galaxies which have lost most of 
their disc and so transformed to SO morphological type, but 
still have only a small bulge which can be cleared easily and 
rapidly by the AGN. 



4.2 Gas replenishment in clusters 

The calculations above are valid for field galaxies, where 
internal processes of gas depletion dominate. If a galaxy falls 



5 DISCUSSION 

We have shown that supermassive black holes residing in the 
bulges of spiral galaxies should have an M — a relation that 
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can be approximated by a power law slightly steeper than 
M oc a 4 , if gas depletion due to star formation is considered. 
Typically, such galaxi es would have M < M CT as predicte d 
in our earlier papers l|Kind 12001 |2010t iKing et al.ll201ll ). 
During an event triggering strong central accretion, such as 
a major merger, the SMBH is likely to grow by a further 
factor of < 7.5 (i.e. for < 2 Salpeter times) before it can 
drive the gas out of the galaxy halo. This factor is an upper 
limit, because we calculated it for an isothermal background 
potential, which is the most difficult to escape from. 

One might think that since the SMBH is not imme- 
diately aware of the processes occurring around the virial 
radius of the halo, it could grow for longer than the calcu- 
lated ~ 2 Salpeter times, as long as the gas reservoir feeding 
it is not depleted. However, the total mass of gas that must 
be accreted in this process is of order 1O 9 M0, perhaps even 
more. This is too high to reside in a single accretion disc, and 
the SMBH is very probably fed in a series of stochastic feed - 
ing events (|King fc Pringlell2007l ; iNavakshin fc Kinell2007h . 
Feedback therefore affects the SMBH feeding reservoir, and 
so it seems more likely that the reservoir is fairly efficiently 
depleted. While gas may fall back on to the SMBH if it has 
not yet escaped the halo and trigger subsequent bursts of 
activity, it cannot do so once the halo is depleted, so the 
AGN should switch off very soon after the outflow clears 
the galaxy. 

Our results suggest that there should be a difference 
between the M — a relations for spiral star-forming galax- 
ies and for red and dead ellipticals. The latter should have 
SMBHs that are systematically overmassive when compared 
with the prediction of the M — a relation for the spiral galax- 
ies. This mass difference is not huge - less than an order 
of magnitude - but it may become apparent when better 
SMBH mass and spheroid velocity dispersion measurements 
become available. 

There is another possible offset in the M — a relation 
depending on the galaxy environment. Galaxies in clusters, 
especially those residing near cluster centres, may accrete 
gas from the hot cluster halo via cooling flows. This pro- 
cess may balance the depletion by star formation at a gas 
fraction ~ 0.2, preserving the old M — a relationship (eq. 
in galaxies that had been spirals by the time they en- 
tered the cluster environment; even though they lose most 
of the gas from the outskirts, the central spheroidal compo- 
nent remains unperturbed and the analysis in Section 13.21 
still applies. Cluster ellipticals would also be a factor < 7.5 
more massive. On average, cluster galaxies have black holes 
slightly more massive (by a factor ~ 2.5) than their field 
counterparts. Our findings are consistent with the recent 
discovery of two SMBHs in the centres of giant elliptical 
galaxies that have masses in excess of 10 10 Mq. They are 
both Brightest Cluster Galaxies and our model predicts that 
their masses should be approximately an order of magnitude 
above the relation derived in earlier works. 

In summary, we have found that there may be three (or 
possibly even four) M — a relations in total, depending on 
galaxy morphology and environment: 

• spiral galaxies with evolved bulges have Mbh ~ 1.3 x 
lO 8 (72(^fM0 with j3 < 1 depending on star formation effi- 
ciency and the relation between the bulge size and a; 

• spiral galaxies with evolved bulges residing in gas-rich 



cluster environments have Mbh ~ 3.7 x lO 8 <rf oo M0 (the 
original M — a relation); such galaxies, however, may be 
extremely rare due to high merger probability in clusters; 

• field elliptical galaxies have Mbh ~ 9.8 x 10 s a 2 oo^-B> 

• elliptical galaxies close to cluster centres have Mbh ~ 
28 x lO 8 cr| oo M . 

There is some uncertainty in the slope, especially for 
field galaxies. It is clear that, as bigger bulges have larger 
velocity dispersions, the slope should be steeper than 4, but 
its precise value depends strongly on the assumed relation 
between Rb and a (both its slope and intercept). 

The differences among the four relations are small, at 
most factors of a few in the normalization. Intrinsic scat- 
ter in each population may blur the distinctions further. In 
addition, since elliptical galaxies tend to have higher veloc- 
ity dispersions than spirals, and cluster galaxies than field 
galaxies, the overall effect of combining the relations may be 
to increase the observed slope of the M — a relation. We illus- 
trate this in Figure [2] We plot the four M — a relations (as- 
suming ft = 0) over a polygon representing the approximate 
locus of data points taken from Figure 3 of lMcConnell et ahl 
(2011). All four lines intersect the locus. The large fraction 
of the locus being below the four lines is reasonable since the 
predicted relations are only upper limits for different envi- 
ronments. Furthermore, there is a hint that higher values 
of a correspond to higher M — a relations. We also plot a 
line with M oc a 6 , which roughly bisects the observed lo- 
cus and so might emerge from a fit to the data unsorted by 
galaxy type, especially in the a > 100 km/s part, where we 
do not expect contamination from nuclear st ellar clusters 
^McLaughlin et alJl20od : lFerrarese et alj|2006ft . 

There have been recent tentative claims that the M — a 
relation may be slightly different for different galaxy types, 
and that th e general relation may show an upturn at highe r 
values of a l|Marconi fc Hunj|2003l : iMcConnell et al.ll201ll l. 
Our considerations here suggest that both features may be 
natural consequences of wind feedback. 
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